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Abstract 

We extend concepts from microwave engineering to thermal interfaces and explore the principles 
of impedance matching in ID. The extension is based on the generalization of acoustic impedance 
to non linear dispersions using the contact broadening matrix r(w), extracted from the phonon 
self energy. For a single junction, we hnd that for coherent and incoherent phonons the optimal 
thermal conductance occurs when the matching r(a;) equals the Geometric Mean (GM) of the 
contact broadenings. This criteria favors the transmission of both low and high frequency phonons 
by requiring that (1) the low frequency acoustic impedance of the junction matches that of the 
two contacts by minimizing the sum of interfacial resistances; and (2) the cut-off frequency is near 
the minimum of the two contacts, thereby reducing the spillage of the states into the tunneling 
regime. For an ultimately scaled single atom/spring junction, the matching criteria transforms to 
the arithmetic mean for mass and the harmonic mean for spring constant. The matching can be 
further improved using a composite graded junction with an exponential varying broadening that 
functions like a broadband antireflection coating. There is however a trade off as the increased 
length of the interface brings in additional intrinsic sources of scattering. 
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I. INTRODUCTION 


Thermal management of devices at the nanoscale requires a comprehensive understanding 
of interfacial phonon transport. At length scales on the order of tens of nanometers, phonons 
flow quasi-ballistically and scattering at interfaces plays a dominant role on the overall 
thermal conductance [l-5|. Many applications stand to benefit from decreasing interfacial 
scattering; for instancy the heat generated in current transistors is blocked by their low 
thermal conductance g, [Tj], raising device temperature and impacting performance and 
reliability js]. One way to increase interfacial conductance is by inserting a thin junction 
layer that allows a gradual change of material properties and bridges phonons across the 


enlargement on 
Nevertheless, 


junction. A recent experimental study on Au/Si interfaces showed a 2 folc_ 
interfacial thermal conductance by the addition of a thin Ti layer (ps 7nm) 
the principles to choose the junction material that maximizes interfacial thermal conductance 
remain unclear. Finding those principles is not an easy task, since they need to achieve the 
largest transmission for a broad band spectrum of phonons, with non-linear dispersion and 
well defined translational and rotational symmetries. 

Some studies have already provided clues on junction properties that correlate with in¬ 
creasing interfacial conductance. It was suggested that maximizing the density of states 


overlap maximizes thermal conductance jd, Il0l-ll2|. Also, it was proposed that maximum 


conductance requires the junction energy flux to equal the Geometric Mean (GM) of the 


contact fluxes 


13|. For atomic junctions, the maximum happens for the Arithmetic Mean 


AM) of the contact masses and the Harmonic Mean (HM) of the contact spring constants 
IJ, [iSj. In spite of the many advances towards hnding the junction material that maximizes 
interfacial thermal conductance, a unifying and quantitative physical picture is still missing. 

In this paper we explain how to choose the mass, spring constant and length (m, /c, L) of 
a single or several junction materials to maximize the average of phonon transmission over 
frequency (T)^, which is equivalent to maximizing the thermal conductance of ID chains 
(Fig. [T]). Our analysis extends the acoustic impedance to the non-linear dispersion limit 
using the contact-induced broadening matrix F, obtained from the imaginary part of the 
self-energy matrix in the Non-Equilibrium Green’s Function (NEGF) formalism jl^ . 

For a single junction material, we hnd that the desirable broadening is the Geometric 
Mean (GM) of the contact broadenings. In the coherent limit, this choice maximizes the 
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FIG. 1. Interfacial thermal conductance is maximized by inserting a graded junction with an 
exponentially varying “impedance” F^ = F/e^” (Eq. [T9l) . However as the length between the 
contacts increases, the conductance gained by “impedance” matching is dominated by losses due 
to other scattering mechanisms. 

resonance peaks of the transmission fnnction (Sec. IIII|) . while in the incoherent limit it 
minimizes the snm of interfacial resistances (Sec. IIVII . To satisfy the GM condition we need 
(1) to choose the acoustic impedance as the GM of the contact impedances to favor low 
frequency phonons and (2) to set the cut-off frequency close to the minimum of the contact 
cut-off frequencies to favor high frequecy phonon and reduce the spillage of states into the 
tunelling regime (Fig. [3h). In Sec. IIIIl those sub conditions are translated into rules to 
choose the junction mass, spring constant, group velocity, Debye temperature or density of 
states. As the junction becomes a single atom or bond, the resonances are pushed outside 
the allowed frequency interval and the GM condition transforms to the Arithmetic Mean 
(AM) of the masses or the Harmonic Mean (HM) of the spring constants of the contacts. 

For multi-junction materials, the GM condition evolves into an exponentially varying 
broadening, since the GM of two quantities stays near the lower partner. That variation 
generates a broadband anti-reflection coating in the coherent limit (Fig. |3h Sec. |V]) and min¬ 
imizes the sum of resistances between dissimilar materials in the incoherent limit (Sec. IVip . 
In this way, thermal conductance is pushed to the maximum physical limit allowed by the 
contacts as the number of junctions increases. Finally, we discuss the trade-off between 
increasing conductance by adding junctions to improve the matching of the contacts and 
decreasing conductance as a result of enlarging the interface and adding intrinsic scattering 
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(Fig. [I]). 

The results of this paper can be used for modeling thermal conductance of Self-Assembled 
Monolayers (SAMsland for extending ID propagation models of long-wavelength phonons 
in semiconductors [l7|, ll8|. It also provides a solid ground to expolore the 3D maximization 
problem since a 3D system can be decomposed into decoupled “ID like” subsystems, with the 
caveat that in each of them the nature of the polarization vectors and their overlap should 
play a role, depending on the symmetry breaking disorder at the interface. The exponentially 
varying broadening provides one example of phonon transmission engineering inspired by 


microwave engineering of broad band hlters 


19| . Similar analogies can be found for the 


binomial or Chebyshev matching transformers, which optimize the flatness or bandwidth of 
the transmission function respectively 19|. We can further envision engineering an ultra 
narrow transmission that hlters particular phonons or an extremely low transmission to 
improve thermoelectric hgure of merit. 


II. METHODOLOGY 


Thermal conductance is dehned as the ratio between heat current and tempera' 
difference across the contacts AT, which according to Landauer formalism is given by 

G® = 


ure 


20 | 


-tujjMT- 


2tt 


(9T’ 


( 1 ) 


where ho; is the energy of a phonon, M is the number of conducting channels, T is the 
average transmission per channel and N the Bose-Einstein distribution. The factor MT can 
be found using NEGF formalism from MT = Trace {riGr 2 Gl}, with G the retarded Green’s 
function of the system of interest and Tj the broadening matrix, which is the anti-hermitian 
part of the self-energy and describes the interaction of the system with contact i j^, 221. 
The number of conducting channels depends on the dimension and force constants of the 
system. In particular, for a ID system with 1st neighbor interactions M reduces to one. We 
further simplify Eq. [1] by assuming a temperature larger than the Debye temperature, so 
the factor hwON/dT is almost constant over the allowed frequency spectrum. In this way, 
G^ oc (T)(^ and maximizing (T)^^ maximizes (Table H]). 

In spite of our second assumption to simplify G^, our results are not limited to high 
temperatures. As we explain in Sec. IIIII and HVl the necessity of Ending a robust maximum 
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requires us to maximize ((T),^)^ the frequency average of the phase average of T. The phase 
average transforms T into a slow varying function by washing away the interference wiggles. 
Therefore, for the maximization process we can pull {T)^^ the phase average of T out of the 
integral and Eq. [T] can be approximated by 

^macc 

« ((r)^)„ I (2) 

0 

Now that the temperature dependence has been isolated, maximizing ((T)^)^^ yields to the 
more robust maximum of G^. The approximation in Eq. [2] breaks down for very low tem¬ 
peratures because the frequency average should be done only over the exited phonons which 
constitutes a frequency range shorter than [0,u)max]- Nevertheless, the material properties 
that maximize ((T)^)^ only differ from the exact solution for temperatures below 10% of 
the Debye temperature. 


T 

Coherent transmission 

(T). 

Frequency average of T over {Q-,00max] 

(T)^ 

Incoherent transmission or phase average of T 


Frequency average of the phase average of T 


TABLE I. Summery of the notation for the different types of phonon transmission used in the 
document. The maximum freuquency Umax is defined by the minimum cut-off frequency of the 
contacts. 


III. SINGLE MATERIAL, COHERENT. 

GEOMETRIC MEAN AND ANTI-REFLECTION COATING. 


Analogous to the quantum mechanical transmission of a particle through a barrier, the 
phonon transmission function T for the system in Fig. [2^ is found from the ratio between 
the incident and transmitted heat currents. After assuming plane waves solutions, imposing 
boundary conditions at the interfaces and replacing the acoustic impedance Z with the 
broadening matrix T 16|, T is given by 

47ir 


T = 


{jir -b 1)2 cos2(giL) -b (7zi -b 7ir)^ sin2(giL) ’ 


( 3 ) 
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where the wave vector of the junction qi comes from the dispersion relation, i.e. uj = 
ujci |sin(gia/2)| with the cut-off f req uency Ud = 2^Jkl/ml, and 70 ,^ = Fa/r^. The broad¬ 
ening function F for a ID chain j23| describes how easily phonons can move from one site 
to the next and is related with the acoustic impedance through [16 1 


F = 2ujpVg{uj) = 2ujZ{uj), 


(4) 


with p the mass density and Vg the group velocity. Eq. [3] is valid for frequencies smaller 
than all the cut-off frequencies. When cuci < OJ < Uci,ujcr, phonons tunnel through the 
junction and the cos and sin become their hyperbolic counterparts. When the frequency is 
larger than the maximum frequency u > m\\i{uJcuUJcr) = ojmax) T = 0 dehning the allowed 
frequency range. 

Note that the transmission function T also follows from adding coherently all the wave 
paths that cross the interface 


T = Jrl \tf =lrl \tlltlr + klTIrT H- f (5) 

_ TiiTir 

1 -I- 2\/RirRii cos(2giL) -|- RirRu ’ 

where at each abrupt interface, the reflection and transmission coefficients are = ( 7 q,^ — 
l)/(7a/3 + l) and tap = (27q,/3)/(7q,/3 + 1 ), and the reffection and transmission are Rap = \Tap\^ 
and Tap = 7/3«IWP- 

As shown in Fig. [2b, T is an oscillating function bounded by the transmission of the 
abrupt interface between contacts 


Tir — 




( 6 ) 


which is prespecified, and an envelope function that depends on our choice of junction 


env — / N r> • 

(7n + 7 ir)^ 

Therefore, as long as Fi lies in between F; and F^, T^nv > Tir, (T)^ > (Tir)^ and the ballistic 
thermal conductance of the gradual interface is larger than the abrupt interface. 

Similar to an anti-reflection coating, we can eliminate reflection at a single frequency by 
destructive interference of the wave path undergoing a single reflection with the sum of paths 
undergoing more than one. The reflection coefficient r is given by the sum of all the wave 
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FIG. 2. a) ID crystal modeled as a chain of masses joined by springs, b) and c) The choice that 
maximizes {T)^ is an interplay between pushing the envelope up and fitting the right number of 
resonances inside the allowed frequency range. The figure shows an example where the largest 
envelope generated by the GM of contact broadenings (b) does not guarantee maximum {T)^ (c). 
It also shows the transmission for the abrupt interface between contacts and the incoherent 
transmission (T),^, which washes away interference (no wiggles). 

paths that return to the incident contact 

r = rii + ^ - 

= ^'1 + -STl • 

To eliminate reflection R = \r\^ = 0, the factors in Eq. [8] must have equal magnitude and 
opposite phase. The hrst condition requires the junction broadening to be the GM of the 
contact broadenings Ti = -y/TiTV, and the second condition happens if the length of the 
junction Li hts exactly a quarter wavelength. 

The choice that maximizes (T)^ is an interplay between pushing T^nv up and htting 
the right number of resonances, since our choice of (mi,fci,Li) affects both the envelope 
function T^nv and the position of the resonances. Intuitively one might think that maximizing 
the envelope over the allowed frequency interval, by making Ti the GM of the contact 
broadenings, generates the largest (T)^. However on a hnite frequency range, the average 
fluctuates according to the number of resonances within the interval, with larger fluctuations 
for fewer resonances. For instance, the maximum (T)^ happens for a single resonance around 
bJmaxj (3/2 — 2/37r^) (Fig. |2fc). This maximum is not very robust since it relies on our ability 
to place a single resonance at a particular frequency, which is limited by our experimental 
accuracy on mi, ki and Li. 
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A more robust condition consists of choosing Fi as the GM of the contact broadenings 
over the allowed frequency interval. This condition maximizes the integral over cu of the 
phase (v9 = qiL) average of the transmission (See Eq. [T3]) 

^GMi’^GMjkcM) max / du\/Tj^\/T^, (9) 

0 

which is equivalent to pushing the envelope function up (Fig. |2]3). Although choosing the 
GM condition does not guarantee the maximum (T)^, its difference with the maximum is 
bounded by — {Tqm)^ < 0.1(1 — {Tgm)^^)- But more importantly, the condition is 

more robust because it does not rely on interference effects. 

The broadening for the GM condition Tqm favors the flow of both low and high frequency 
phonons by making Fi close to vTiTV over the entire frequency range (Fig. [Sb). Note that 
the GM of acoustic impedance by itself does not guarantee high (T)^ (Fig. EH)- Indeed, if 
our choice of mi, ki makes cUd < oomax, phonon transmission is drastically cut by tunneling. 




FIG. 3. (a) Choosing the impedance Zi = y/ZiZf does not guarantee high (T)^ because if the 
cut-off frequency lies below ujmax, tunneling cuts down transmission (shaded region). On the other 
hand, choosing the broadening Fi Ri yTTk (b) fixes a unique pair mcM^kcM that favors the 
flow of both low and high frequency phonons. The figure shows a similar result for the incoherent 
transmission average {{T)^)^]. Note that F for the maximum (T)^ {flmax) intersects -v/TTV where 
the single resonance is located (Fig. [IJi). 

A back of the envelope estimation for m^M and kcM follows by graphically overlapping 
Fi(mi,fci) with a/T/TV (Fig. Eb). For low frequencies the slope of F is dictated by the 
impedance, so to favor the flow of low frequency phonons we want Zi = \fZf7fl. To favor the 
flow of high frequency phonons, we equate the frequency cu* at which a/T/TV is maximum to 













the frequency at which T is maximum cjci/\/2 = oo*. Then, from Z = y/mk and ojc = ‘2y/k/m 
we solve for uigm and kcM 

iij(^\Z\ ‘2.Z\ 

kcM ~ —— mcM ~-• ( 10 ) 

Z Ucl 

From Eq. [10] we can also infer how to choose the junction group velocity Vg = ay/k/m and 
the Debye temperature To = fyuJo/kB with ojd = oJcT^/2 and ks the Boltzmann constant. 
Using Eq. H] we can also express the GM condition in terms of the density of states per unit 
cell {DOS) 

( 11 ) 

irDOS’ ^ ’ 

which quantihes how to choose the overlap of density of states to maximize thermal conduc¬ 
tance. 

When the junction material becomes a single mass or bond, the GM condition to maxi¬ 
mize (T)^ transits to the AM of the contact masses or the HM of the contact spring constants 
respectively. This transition happens when ki increases enough or mi decreases enough to 
push Tenv below Tir (Fig. [2|) making (T)^ < Then it is better to have an abrupt 

interface between the contacts which take us back to the AM or HM condition. In par¬ 
ticular, if the spring constant ki » uj‘^^^mi{Lja)"^, then sin^(giL) /kia and 

Fi 2ujy/kimi. Replacing the approximations into Em 3] and letting ki tend to inhnity we 
recover the transmission for a single atomic junction 16] 

1 ^ 

The maximum of Eq. [12] happens when mi = 0, which recovers the AM condition for the 
atomic mass linking the contacts (Fig. [2^). Another way to rationalize the transition from 
the GM condition is by realizing that increasing ki or decreasing mi pushes the resonances 
out of the allowed frequency interval. Then, without resonances available to reach T^nv, 
maximizing T^nv with the GM condition is no longer useful. 


IV. SINGLE MATERIAL, INCOHERENT. 
STILL GEOMETRIC MEAN. 


When the phase gained by a phonon in its transit between interfaces is randomized, 
interference disappears and we can think of phonons as classical particles. The phase {ip = 
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qiL) average of the transmission in Eq. [3] or the incoherent transmission is given by 
(Fig. ED 


(F) = 


47«r 


= 


(13) 


(7;r- + l)(7n + 7ir) 

Since Ti^ is prespecihed, this expression is once again maximized by the GM condition, for 
which Tenv ~ 1 and {TGM)ip ~ \/Tw > Tir (Fig. Eb). Note that Eq. [T3]is the same expression 
obtained by adding incoherently all the possible wave paths from the left to the right contact 


(F) — TiiTir + TiiRirRilTir 


TnTi 


11 ^ Ir 


(14) 


1 RllRlr 

with Tajs = ^'yajs/ilajs + 1)^; fFe transmission at the interface between material a and /3, 
and Rajs = 1 — Tap the reflection. 

The GM condition can be rationalized as a minimization of the sum of two opposite 
resistances, which are proportional to the ratio between phase average reflection and phase 
average transmission 2^. Reorganizing Eq. [14] to expose the additive property of resistance 

(F), 1 - (F)^ 


25| 


(F), 


1 — F^i ^ 1 — Tir 


(15) 


(F)^ Tn ’ 

it follows that a variation of the junction produces opposite trends on the single interface 
resistances. Therefore, Eq. [T5] is minimum when the individual resistances are equal, which 
leads to the GM condition. 

Phase randomization can arise from lack of experimental control. For instance, the mea¬ 
sured thermal conductance on a SAM junction is the sum of transmissions over molecules 
with slightly different lengths determined by the cross-sectional variation of the SAM. Ac¬ 
cording to the variation strength, the average transmission can be approximated by 

sin(2(5) 




25 


+ (F)., 


(16) 


with 25 a uniformly distributed phase interval around over which we average T. This 
expression follows from finding the effect of phase randomization on the peak of a sinusoidal 
function and using it to interpolate the two known limits. Note that 5 = 0 recovers the fully 
coherent limit Eq. [3] while 5 = tt/2 the fully phase randomized limit Eq. [131 

In contrast to the coherent limit, the incoherent limit does not provide a condition for 
the junction length. To be consistent with our assumption that the phase is randomized 
between interfaces, we need to guarantee that Li is longer than the phonon phase coherent 
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length. In case the phase randomization happens at the interfaces, Li is limited by the 
longest distance between interacting unit cells. That is the smallest block of material that 
allows us to define impedance. 


V. N MATERIALS, COHERENT. 
EXPONENTIAL. 


For several junction materials in between the contacts (Fig. |3^), enforcing the GM con¬ 
dition for every three consecutive pieces translates into an exponentially varying broadening 
(Eq. [19]). This trend follows from taylor expanding F(a:) = to sec¬ 

ond order and solving the corresponding differential equation (F')^ = FF". As we increase 
the number of materials, this variation generates a banded transmission that monotonically 
widens. In this way, we approach the maximum limit of thermal conductance imposed by 
the contacts. The idea was inspired by the exponentially tapered impedance coupling used 
to design broad band transmission lines [19|, which constitutes just one example of the 
extrapolation possibilities from Microwave engineering to phonon engineering. 

The transmission for a system with N junction materials is a complicated function dic¬ 
tated by interference patterns. However, to maximize transmission, we intuitively expect a 
small reflection at each interface so we can approximate the reflection coefficient r by the 
sum of wave paths undergoing a single reflection, i.e. 


r K.rii+ -f- r23e*^‘^^ -f- 


+ rjvre' 




( 17 ) 


At this point we have N exponential basis vectors to £t a desired reflection function, which 
is a rich problem with different answers according to the design criteria for the function. 
We further simplify the expression by using our intuition of anti-reflection coating, where 
to eliminate the first reflection rn we need r — rn to have the same magnitude and opposite 
phase. Then we assume equal phase gained between interfaces, i.e. ipn = 2nip, to get complex 
roots of 1 equally spaced in Eq. [T3 and equal reflection at each interface rj_j+i = rj+i j _|_2 = 
to assure equal magnitude. Then Eq. [T7I turns to 


N 


•E- 

j=0 


r ^ Vi > = Ti 


]_ _ g*2(7V+l)(/j 
2 _ ^i2(f 


( 18 ) 
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As N increases, the resonances (r = 0 -H- = 1), the anti-resonances = —i) 

and the decreasing trend of generate a reflection that broadens and tends to zero or 
equivalently a transmission that broadens and tends to unity. Fig. 0)3 shows the widening 
of the transmission when we increase the number of junctions from iV = 1 to iV = 5. 


^1 Left Contact Junction Materials Right Contact 
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FIG. 4. a) ID crystal modeled as a chain of masses joined by springs, b) Transmission T and its 
phase average (T)^ for contacts coupled through an exponential variation of broadenings with a 
single junction N = 1 and with N = 5 junction materials. As we increase the number of materials 
transmission widens increasing (T)^. The spatial variation of acoustic impedance is shown in c) 
and of cut-off frequency in d) for N = 1 and N = 5. The crosses show the back of the envelope 
calculation given by Eq. [20l 


Equal reflection at each interface follows from imposing the GM condition on every three 
consecutive materials, which leads to an exponential variation 

N-\-l — n n 

F„ = F, (19) 

with ( = ln(F;/Fr)/(iV -|- 1). This variation maximizes the frequency average of the inco¬ 
herent transmission (See Eq. [23] and below). Therefore we hnd the and kn that better 
satisfy Eq. [13] by maximizing ((T)^)^. Similar to Sec. IIIIl a useful back of the envelope 
approximation follows by graphically overlapping the desired and possible broadenings. To 
match the low frequency slope we need Zn = a^d to match the high 
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frequency spectrum we need Ucn = with cu* the frequency at the maximum value of 

Eq. [191 Then and kn are approximated by 




•'cn^n 

^- nrir, ~ 




00c. 


( 20 ) 


and rules for the variation of the group velocity, Debye temperature and overlap of density 
of states can be found following Sec. IIIII Fig. litjd show the spatial variation of junction 
impedances and cut-off frequencies for = 1 and N = 5. Note that the trend of impedance 
follows an exponential variation that favors flow of low frequency phonons while the trend of 
cut-ojf frequencies is stable around oomax to favor flow of high frequency phonons. 

Equal phase gain between consecutive interfaces needs 

qnLn = —2sin”^ (21) 

® V 00 cn J 

which we can only satisfy at a single frequency uo**. A careless choice of oo** and ip might lead 
to impractically small lengths (L„ ~ 2a) or long lengths where phonons are not coherent 
{Ln ~ mean free path). A sensible choice is for example uo** = uOmax/‘^ and = 27r, which 
generates lengths on the order of ~ 10a ~ 5nm. 

The theory of small reflections shown in this section is valid as long as we can neglect 

the contributions from paths that involve more that one reflection, i.e. |rj|^ << |rj|. For 

Ti < Tr and if we consider |rj| <0.1 small enough, the minimum number of junctions A^m 
needed to be in this regime is 

- 1 . ( 22 ) 

For instance, for a 10-fold mismatch r^/T/ = 10 Nmin ~ 2.5 and for a 100-fold mismatch 


VI. N MATERIALS, INCOHERENT. 

BACK TO EXPONENTIAL. 

When phonon phase is randomized in between interfaces, the exponentially varying broad¬ 
ening minimizes resistance and generates maximum Total resistance decreases as 

the number of junction materials N increases because the addition of resistance coming from 
new interfaces is dominated by the decrease on interfacial resistance due to less mismatch. 
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Nevertheless, as we increase N, the length between contacts increases and other scattering 
mechanism become important. Therefore, there is a sweet spot for N at which we stop 
increasing thermal conductance by adding junction materials. 

Transmission when phonon phase is randomized or incoherent transmission (T)<^ can be 
found from the total resistance of the system similar to Eq. (TS] 22, l2J, [25 1 

1 - (T)^ 


(T), 


1 ~ Tj 

2-^ T- 


(23) 


Eq. [23] becomes minimum when all the factors on the right hand side are equal. This condi¬ 
tion needs the ratios of broadenings between consecutive materials to be the same and leads 
to the exponential variation in Eq. [19] (Fig. [Dd). Replacing that variation into Eq. [23] give 
us the maximum (T)^, which can be approximated by 

4iV 


(T) 


V’ i„2 


ln^(R/r,)+4iV 


when N tends to inhnite and by 


(T),« 


(24) 


(25) 


for small N. As shown in Fig. [5^, (T)^ increases as N increases. This non intuitive result 
happens because the increment in resistance by the addition of an interface is dominated by 
the decrease of the interfacial resistance of all the interfaces due to less material mismatch. 



FIG. 5. a) Using the exponential variation of broadenings, as the number of junction materials 
N increases, {T)^ and {{T)^)i^ approach to one. b) As N increases the length between contacts 
increases and phonon flow is affected by scattering mechanism different from interfacial scattering. 
Therefore, there is a sweet spot for N at which we stop increasing thermal conductance by adding 
junction materials. 

As the number of junction materials increases, the length between contacts also increases 
and other scattering mechanisms besides interfacial scattering become important. We can 
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combine all of them using Matthiessen’s rule to obtain an effective mean free path (mfp) in 
terms of the mfp for interfaces Aj and the mfp for other scattering mechanisms As 


A, 


eff 


1 

A“ 


and an effective transmission 


{TK" = 


A 


1 

eff 


L + Xeff 

Note that Eq. [26] and [23 can be derived from equation [23] 


(26) 


(27) 


22 


25[ |. We dehne the interfacial 


mfp for N materials following the exponential variation A* by replacing Eq. [THl into Eq 
dehning the distance between interfaces as Li = L/N and comparing with Eq. [57] so that 


{T)^ = 




L + Aj 


and 


A* = 


LLi 


L -\- Lj 


Ri 


(28) 


Aj is dominated by L, since we expect Li < L and is inversely proportional to the resistance 
at each interface. Thus, as N increases Aj also increases. However, the effective mfp is 
bounded by As. Thus, there is a sweet spot for the number of junction materials at which 
adding more interfaces does not increase the effective mfp and {T)^R is dominated by the 
increment of L (Fig. Od). 


VII. CONCLUSION 

We show how to choose the properties of a single or multiple junction materials between 
prespecified contacts to maximize thermal conductance in ID systems within the coherent 
and incoherent regimes. For a single junction our choice should be the GM of the contact 
broadenings, which requires two conditions: (1) the impedance equal to the GM of the 
contact impedances to favor low frequency phonons and (2) the cut-off frequency close to 
the minimum of the contact cut-off frequencies to favor high frequency phonons. Those 
sub conditions translate into rules to choose the junction mass, spring constant, group 
velocity, Debye temperature or density of states. These rules allow us to pinpoint the role 
of each quantity in bridging phonon across the interface. For multi-junctions materials we 
show that the GM condition evolves into an exponential variation of broadenings. This 
variation pushes thermal conductance to the maximum limit allowed by the contacts as the 
number of materials increases. We also integrate interfacial scattering with other scattering 
mechanisms. We hnd a sweet spot for the number of material caused by the interplay 
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between smoothing the interface by impedance matching and enlarging the distance between 
contacts. 
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